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The microscopic spectral density of the Hermitian Wilson-Dirac operator is computed numerically 
in quenched lattice QCD. We demonstrate that the results given for fixed index of the Wilson- 
Dirac operator can be matched by the predictions from Wilson chiral perturbation theory. We test 
successfully the finite volume and the mass scaling predicted by Wilson chiral perturbation theory 
at fixed lattice spacing. 



I. INTRODUCTION 



For a long time, the Wilson construction of Dirac 
fermions on the lattice was somewhat overshadovifed by 
other formulations with either explicit chiral symmetry 
or remnants thereof. Recently, with the help of improved 
simulation techniques and more powerful computers, nu- 
merical simulations with Wilson fermions have picked up 
momentum again. The initial concern about the explicit 
breaking of chiral symmetry in this formulation can be 
dealt with in precisely the way originally envisaged: by 
simply pushing simulations into the regime of very light 
fermions, even to the physical point of almost massless u 
and d quarks. This requires tight control on the small- 
est eigenvalues of the associated Wilson-Dirac operator 
Dw- If the lattice spacing is not small enough, the 
gap associated with the Hermitian Wilson-Dirac operator 
D5 = 75(DvF + 'm), where m is the quark mass in con- 
tinuum language, shrinks due to small eigenvalues pop- 
ulating the region between —m and m. This can cause 
severe numerical instabilities in the simulations. Studies 
of the spectral gap and its potential disappearance have 
indicated that simulations may be continued to the phys- 
ical point even at realistic lattice spacings of present-day 
computational power [l[ . 

Sparked by this renewed interest in simulations with 
Wilson fermions on the lattice, also work on analyti- 
cal approaches to understanding the chiral limit of these 
fermions has intensified. As is always the case, it is ad- 
vantageous to phrase the discussion in terms of the perti- 
nent low-energy effective field theory, Wilson chiral per- 
turbation theory (WCPT) [2|. Within this framework 
Sharpe [| was the first to consider the effects of finite lat- 
tice spacing, a, on the spectrum of the Hermitian Wilson- 
Dirac operator. More recently, a complete description 
of the microscopic details of the Wilson-Dirac operator 
spectrum up to and includiiig corrections has been 
presented for both quenched [l-Q and unquenched cases 



[1, Q . Here we will test these recent predictions against 
quenched lattice data. 

The microscopic description works, at leading order in 
the chiral counting rule of the finite-volume e-regime, at 
two equivalent levels: (i) a Wilson chiral random matrix 
theory (WRMT) and (ii) the zero-mode sector of WCPT 
to order a^. Beyond leading order only the approach 
based on WCPT survives. 

The 75-Hermiticity of the Wilson-Dirac operator. 
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(1) 



together with the explicit form of Dw tells us that the 
eigenvalues of Dw come in complex conjugate pairs jioj . 
a feeble remnant of chiral symmetry. In addition, 75- 
Hermiticity of Dw ensures that the Hermitian Wilson- 
Dirac operator. 



D5 = IbiDw + m), 



(2) 



as the name suggests, is Hermitian. While the eigenval- 
ues of Dw can be computed by, say, the Arnoldi algo- 
rithm, it is clearly advantageous to focus instead on the 
Hermitian Wilson-Dirac operator D5 , for which fast algo- 
rithms exist that provide an ascending sequence of lowest 
eigenvalues. 

Here we test the predictions for the quenched micro- 
scopic spectral density of D5 that were presented in 0, 
against lattice data. All results of [1,[5[ are given for fixed 
index of the Wilson-Dirac operator. This index v is de- 
fined for a given gauge field configuration by 



^ sign((fc|75|A:)). 



(3) 



Here, |fc) denotes the fc'th eigenstate of the Wilson-Dirac 
operator. Only those eigenvectors which correspond to 
a real eigenvalue of the Wilson-Dirac operator have a 
nonvanishing expectation value of 75 [lo| and hence con- 
tribute to the index. 
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In the microscopic Umit the real eigenvalues of the 
Wilson-Dirac operator are predicted to be located in re- 
gions with width of order 1/F H S, 0, where V is the 
four volume, and we choose only to include the physi- 
cal/leftmost of these regions in the above definition of 
the index. This is indicated by the prime on the sum in 
Eq. 

The predictions Q for the quenched microscopic 
spectral density of are based on the e-regime of Wil- 
son chiral perturbation theory. At leading order in the 
e-regime of WCPT three additional low energy constants, 
We, W7 and Wg parametrize the leading corrections due 
to the nonzero lattice spacing (order corrections) 0- 
There are large- iVc arguments {Nc is the number of col- 
ors) to the effect that the Wilson low-energy constants 
Wg and W7 may be small [1] . For the case Wq = Wt = 
it was argued in 0, Q that only Ws > correctly de- 
scribes lattice QCD with a nonnormal and 75-Hcrmitian 
Dirac operator. In agreement with this wc will demon- 
strate here that the quenched microscopic eigenvalue den- 
sity of Z?5 obtained on the lattice can be matched by the 
prediction from WCPT with the low energy parameters 
We = 0, VK7 = and Wg > 0. 

Wc identify the effects of nonzero Wq and Wr in the 
spectrum of D5. and conclude that they do not change 
the above conclusion. 

Finally, we explicitly compute the real eigenvalues of 
Dw and compare to the predictions. 

Wilson chiral perturbation theory has also been consid- 
ered in the finite- volume scaling limit where the a^-terms 
can be pulled down from the action and treated as small 
perturbations The analytic constraint on Wg that 

wc discussed above is not obvious in such a scaling limit. 

Our paper is organized as follows. In the section just 
below, we first briefly recall those results from WCPT in 
the e-regime that are most relevant for the present lattice 
study. Section IIIII gives the definition of Wilson chiral 
random matrix theory. In section IIVI we give details of 
the lattice simulations we have performed. We compare 
our lattice data to the predictions from Wilson chiral 
perturbation theory in section|Vl Section lVTl contains our 
conclusions. 



II. PREDICTIONS FOR THE MICROSCOPIC 
SPECTRUM 



Here we recall the predictions for the spectral density 
of ^5, 

P5(A^m;a) = ^^^(Ai-A5)y (4) 

The eigenvalues A^, defined by 

D5^k = AfVfc (5) 



are considered in the microscopic limit where V ^ 00 
while 

mV, XlV, aV (6) 

are kept fixed. 

The microscopic spectral density follows from WCPT. 
To 0{a^) accuracy this effective theory was written down 
in expanding in all operators that contribute to this 
order in the Szymanzik scheme. To order the zero- 
mode partition function at fixed index v can be obtained 
Q by decomposing the momentum zeromodc partition 
function according to 

00 

ZNf{m,0;a)= ^ e'^^'^Z]^^ (m; a). (7) 
If— —00 

In this e-regimc, the partition function reduces to a uni- 
tary matrix integral 

Z'^ (m,z;a)^[ dU defU e'^^^l (8) 

Ju{Nf) 

where the action S[U] for degenerate quark masses is 
given by 

S = yEVTr(f/-f C/t) + ^SFTr(C/-C/'f) 

-aVWe [Tr {U + C/^)]^ - a^VWvlTT {U - C/t)]^ 
-a^VWs.TviU'^ + U^'^). (9) 

In addition to the chiral condensate, E, the action con- 
tains the low-energy Wilson constants We, W7 and Wg 
as unknown parameters [13]. In the continuum the de- 
composition (17]) into a sum over sectors of index v is 
unambiguous hj. At nonzero lattice spacing the terms 
of order in the chiral Lagrangian, however, allow one 
to make different decompositions corresponding to dif- 
ferent definitions of 9 and j/ at a 7^ 0. As was shown 
in Ref. Q, the specific decomposition JTHH]) corresponds 
exactly to the definition of v given in Eq. Moreover, 
in the microscopic limit the real modes arc predicted to 
be in well seprated regions each width width of order 
a^Wi/Yj ~ ^/V. The gauge field configurations which 
contribute to the individual terms in the sum ([7]) can 
thus be identified on the basis of the spectral fiow jT3| . 
This provides a firm foundation for choosing this partic- 
ular lattice definition of the index. 

To simplify our notation, we absorb the factor of VWi 
into af and the factor VYi into to, and 

a^VWi a\ raVH rh 

zVE^z A^yE^A^ (10) 

The explicit expression for the quenched microscopic 
spectral density can be derived following the procedure 
described in Ref. [3], now extended to the case of a 
source z for the pseduoscalar density. In addition, much 
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like at nonzero chemical potential [T^ [T^ , one must care- 
fully take into account the convergence of the noncom- 
pact integrals. The resulting analytical formula was ob- 
tained in [1, 



P5(A^,TO;a^) 



1 



where 



Im G'''(-A^m;ai), 



(11) 



/ C f]f) 
ds J —icos{e)e^''^-'> 

X cxp[— 777. sin(0) — S7isinh(s) -t- cos{0) ~ii^z~ ie)cosh(s) 
-h4ag(-i sin(6') + sinh(s))^ + 4a7(cos(6') - cosh(s))^ 



-2a|(cos(26l) - cosh(2s))] 
— sm(t^) + I— smh(s) - 
-4(ai + a?)(sin2(0)+sinh2(.s)) 
-f2a8(cos(20) + cosh(2.s) 4^ e 



777 777 Z Z 

— sin(6') + i— sinh(s) + i- cos(6') + i— cosh(s) 



ie+s 



')+2)-(12) 



Existence of the integrals leads to the constraint Wr> 
(for Wq = Wj = 0) as described in detail in Refs. 0, 
The sign of Wg is also crucial for understanding if one 
approaches the so-called Aoki phase of Wilson fermions 
[l7| . A positive sign of (and Wq = = 0) indicates 
that the Aoki phase will appear as the gap in the spec- 
trum of i?5 closes. The analytic predictions of WCPT for 
the spectral density I?5 allow us to follow in detail how 
the gap disappears as m is taken to zero at fixed lattice 
spacing a. 



III. WILSON CHIRAL RANDOM MATRIX 
THEORY 

As we will make explicit use of it below, we give also 
here the corresponding expressions in the equivalent Wil- 
son chiral random matrix theory formulation 0, Q . The 
idea is to introduce the most general large- A'^ matrix that 
has conjugation properties similar to 75-Hermiticity. Let 
this matrix be 



aA iW 
iW"! aB 



where 



A = A^ and = B 



(13) 



(14) 



are (n + v) x (n + v) and n x n complex matrices, re- 
spectively, and W a rectangular complex matrix of size 
{n + ly) X 77. As in Q we use tildes to indicate quantities 
in Wilson chiral random matrix theory which arc ana- 
logues of those in the field theory. We consider the limit 
N = 2n + i' ^ 00 with miV, zN and a^A fixed. 
The matrix D\y is 75-Hermitian with respect to 



where the first (n + v) entries arc -1-1. and the remaining 
n entries are -1. 

Results based on the Wilson chiral random matrix the- 
ory are expected to be universal, as in the case of the 
usual chiral random matrix theory [l8| . For simplicity 
one can therefore take the matrix elements to be dis- 
tributed with Gaussian weight 



P{A,B,W) = e 



(16) 



The partition function of the Wilson chiral random ma- 
trix theory is thus defined by 



Z'^^. = dAdBdwYl det{Dw + mf + Zf%) P. 
f=i 



(17) 



The matrix integrals are over the complex Haar measure. 

In the limit N 00 the partition function (jl7p coin- 
cides with that of ^ when Wq = Wj = 0. It is straight- 
forward to extend the Wilson chiral random matrix the- 
ory so that also double-trace terms corresponding to Wq 
and Wr are included The correspondence between 

parameters of Wilson chiral random matrix theory and 
Wilson chiral perturbation theory is 



Nrh 



mTiV 



N2 



zYM Na 



a^WsV. (18) 



75 = diag(l,--- ,1,-1,' 



-1) 



(15) 



Once this identification is made, one can equally well use 
the Wilson chiral random matrix theory to derive spec- 
tral distributions. This has been demonstrated explicitly 
in @, Q- Moreover, since no compact expressions are 
presently known for individual eigenvalue distributions, 
it is straightforward to generate such eigenvalue distribu- 
tions numerically using the Gaussian integrals of Wilson 
chiral random matrix theory. This will be used below. 



IV. THE LATTICE SETUP 

Here we give a few details of the lattice simulations. 
For this initial study we work in the quenched approx- 
imation. It is advantageous to use a pure gauge action 
which gives configurations with a fairly unique topolog- 
ical charge, i.e., a pure gauge action that suppresses so- 
called dislocations. Experience from simulations with do- 
main wall fermions, for which the dislocations tend to 
increase the residual chiral symmetry breaking at finite 
fifth dimensional extent L^, suggests that the Iwasaki 
gauge action (2^ would be a good choice. This action is 
fairly well studied and appears to have a rather smooth 
continuum limit. Our lattice parameters were chosen 
using rp/g values from interpolation formulae given in 
[21I |2^ and using ro = 0.5 fm to set the physical scale. 

We used a mixture of heatbath and overrelaxation up- 
dates with, as advocated in Ref. [13, Nor ~ 1.5ro/a over- 
relaxation sweeps for each heatbath sweep. We saved 
configurations separated by 100 heatbath sweeps for our 
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Ensemble size ro/o a [fm] L [fm] # cfgs 

A 2.635 16* 5.37 0.093 1.5 6500 

B 2.635 20* 5.37 0.093 1.9 3000 

C 2.79 20* 6.70 0.075 1.5 6000 



TABLE I: Parameters of the pure gauge configurations con- 
sidered. They are generated with tlie Iwasaki gauge action. 



measurements. We found no detectable autocorrelations 
in the ensembles generated. Experience indicates that 
for studies of the e-regime of chiral perturbation theory, 
physical lattice sizes of i = 1.5 to 2 fm arc most use- 
ful. The lattice size is large enough so that several Dirac 
eigenvalues arc governed by the e-rcgimc of CPT, but 
it is small enough that sufficient statistics in sectors of 
given, small topological charge can be obtained. Some 
relevant parameters of the three ensemble considered in 
this study arc summarized in Table HI 



A. Spectral flow 

In order for the match to Wilson chiral perturbation 
theory to be valid the assignment of index ly to sl given 
gauge field configuration must be obtained from Eq. ([3]) 
using either a) the (small) real eigenvalues of the Wilson- 
Dirac operator, together with the chiralities of the cor- 
responding eigenvectors or equivalently by counting 
the number of net crossings of zero of the eigenvalues of 
the Hermitian Wilson-Dirac operator as a function of the 
mass parameter toq [l3j and weighting their contribution 
by the sign of the slope at the crossing. Since we did 
not have an implementation of the Arnold! algorithm to 
compute the (complex) eigenvalues of the Wilson-Dirac 
operator, including the relevant real eigenvalues, we used 
the spectral flow strategy. 

The index determined from the real eigenvalues of the 
Wilson-Dirac operator, or equivalently from the spectral 
flow, can depend on the range of the real eigenvalues con- 
sidered, or equivalently, the value of —rricut at which the 
spectral flow is terminated. With —rric being the critical 
mass, at which the pion formed from Wilson quarks be- 
comes massless, rricut should be sufficiently larger than 
nic, but smaller than the value at which the doubler 
fermions become light. For our computations we use one 
HYP smearing of the gauge fields [2j] in the construction 
of the Wilson-Dirac operator to help further suppress lat- 
tice artifacts. Then the choice rricut = 1.0 is safely away 
from both nic and the doubler region. With the Iwasaki 
gauge action and its suppression of dislocations, crossings 
in the neighborhood of rricut are rare, but they do occur 
occasionally. Hence a change in rricut would change the 
index on a few configurations. 



Ensemble 


size 




mo 


# configs 
= 0,1,-1 


A 


16* 


2.635 


-0.216 


1246, 1088, 1045 


B 


20* 


2.635 


-0.216 


379, 319, 322 


C 


20* 


2.79 


-0.178 


1172, 990, 988 


C 


20* 


2.79 


-0.184 


1172, 990, 988 



TABLE II: 



V. THE SPECTRUM OF D5: LATTICE 
SIMULATIONS 

In this section wc compare the microscopic spectral 
density of 1)5 as predicted by WCPT for the quenched 
case to the lattice simulations described above. 

The values of m, and a given in the figures refer to 
TO, and a%. For all plots except Figure |8] we have set 
0-6 = 0-7 = 0. 

A. Volume scaling 

Our first test of the analytic predictions for the micro- 
scopic spectrum of Z?5 concerns the finite volume scaling 
at fixed lattice spacing. That is, the scaling of to, A^ 
and ds with V for fixed a. To this end wc consider the 
two sets of data with /3/u, = 2.635 and toq = —0.216, see 
Table [TTl Because of the larger statistics obtained we first 
match the 16^ lattice data with v = Q to the prediction 
(fT^ . The result is displayed in the top panel of Figure 
[1] We observe that it is possible to make a reasonable fit 
to the data. The range over which the predictions match 
the data is limited to the lowest two eigenvalues, beyond 
which the increasing slope of the data indicates that the 
p-rcgime effects set in. Given the best values found in 
the fit we then scale m by 20'*/ 16* and as by 20^/16^ 
to obtain the analytic prediction for the spectral density 
obtained on the 20^ lattice at the same value of the cou- 
pling and lattice mass. This parameter free prediction is 
plotted against the lattice data in the lower panel of Fig- 
ure [TJ We observe that this finite- volume scaling of the 
analytic prediction is consistent with the 20^ data. The 
limited statistics docs not allow us to resolve the detailed 
structure of the analytic prediction. 

Likewise for = 1 we can test the analytic prediction 
with the parameters fixed from the fit to the 16* and 
V = Q data. The result is displayed in Figure [H Note 
that the analytic prediction for — i/ is obtained from the 
one for v by flipping the sign of A^. To make the most 
of our statistics we therefore overlay the data obtained 
for V = 1 and v = —1 flipping the sign of the eigenval- 
ues of the latter. This is done throughout this paper. 
We observe that the general structure of the 16^ v ^ 1 
data is reproduced by the analytic prediction. Also the 
overall features of the 20* v = \ data are followed by 
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FIG. 1: The eigenvalue density of D5 in the sector with v — Q 
for V = 16'' (top) and V = 20* (bottom) both with = 
2.635 and mo = —0.216. The x-axis has been rescaled by EV, 
where E16* — 157.5. In the top plot the smooth red curve 
displays a fit of the analytic prediction for the microscopic 
density of D5 to the 16* data. The fit values are m = 3.5 
and ag = 0.35 (ae = a? = 0). In the plot below the smooth 
red line shows the prediction for the 20* data (m = 8.5 and 
a% = 0.55) obtained from finite volume scaling of the 16* 
values. We observe that the finite volume scaling at fixed 
lattice spacing is well respected. 



the parameter free analytic prediction which in this case 
follows from finite volume scaling. 

We conclude that Wg = W-j = and > is con- 
sistent with the lattice data at fiiw = 2.635 and that 
the volume scaling at fixed lattice spacing predicted by 
Wilson chiral perturbation works well. 



B. Mass scaling 

Our second test of the analytic predictions for the spec- 
trum of D5 concerns the scaling with the quark mass at 
fixed lattice spacing. For this we choose a 20* lattice 
with = 2.79 and two masses mo = —0.178 respec- 
tively Too = -0.184. 



FIG. 2: The eigenvalue density of Ds in the sector with u = 1 
for V = 16* (top) and V = 20* (bottom) both with I3i^ = 
2.635 and mo = —0.216. The x-axis has been rescaled by SV, 
where E16* = 157.5. The smooth red curves displays the 
analytic prediction for the microscopic density of D5. The 
values, 771 = 3.5 and as = 0.35 for 16* and m = 8.5 and 
as = 0.55 are obtained from the fit in the top panel of Figure 
[T]and by finite volume scaling from 16* to 20* {ag = 07 = 0). 



We first consider the toq = —0.184 data with 1/ = 
0. We make a fit and find the best values of to = 3.5, 
Y,V — 216 and ag = 0.35. (The fact that rh and a takes 
the same value as for the 16*^ lattice at fSj^ = 2.635 is 
accidental.) For these values the analytic curve matches 
the data well over the lowest four eigenvalues, see the top 
panel of Figure [31 

Given these values of to and ag we then test the pre- 
dicted density for 1/ = 1 against the data. As can be 
observed from the lower panel of Figure |3] this again 
works well for the lowest four eigenvalues. Finally, we 
consider also the higher mass data at the same value of 
the coupling. Since the lattice mass has changed from 
Too = -0.184 to Too = -0.178 we have Sm = 0.006 
so that SmJ^V = 0.006 * 216 = 1.3. Hence, we predict 
that the 13 =2.79 data with too = —0.178 should be 
matched with to = 4.8 and a = 0.35. As can be seen 
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FIG. 3: The eigenvalue density of Dr, for Pj^ = 2.79, V = 20" 
and mo = —0.184 top: v — and bottom: |!/| = 1. For the 
v = data we find the best values m — 3.5, as = 0.35 and 
"EV = 216. With these values we then test the prediction 
against the |!/| = 1 data. 

from Figure HI this parameter free prediction works well 
for the = as well as for the = 1 sector. Figure 
[5] compares the corresponding parameter free prediction 
for —2 against the data. 

We conclude that Wg = Wr = and > is con- 
sistent with the 20'' lattice data with /3i^ = 2.79 and 
that the mass scaling at fixed lattice spacing implied by 
WCPT works well. For the influence of Wq and M^y see 
below. 



C. The rh and ag dependence 
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FIG. 4: The eigenvalue density of D5 for jSi^ = 2.79, V = 20" 
and mo = —0.178 top: v = Q and bottom: \v\ = 1. We use 
the value = 216 (found for mo = —0.184) and Sm = 
—0.006 to determine the value m = 4.8 for the Pi^^ — 2.79, 
V — 20", mo — —0.178 data. Since the coupling is unchanged 
also ag = 0.35 is fixed by the previous fit. The parameter 
free predictions for u — and |z/| = 1 are then plotted against 
the data. We observe that the mass scaling at fixed lattice 
spacing works well. 



as- Here we demonstrate the strong sensitivity of the 
analytic predictions for the microscopic eigenvalue den- 
sity of D5 with fixed index to the values of m and ag. 
Figure [S] give examples of this strong dependence which 
allow for a precise determination of the best values for 
the match to the data. Figure [7] compares the best fit to 
the continuum, a = 0, curve. In this way we estimate the 
error on ag in our data sets to be at most 0.05 and the 
error on rh to be below 0.5. With increased statistics on 
slightly larger lattices at slightly smaller lattice spacing 
it should be possible to determine m and ag with high 
accuracy. 



Above we have obtained best fits of the analytic pre- 
diction, Eq. (IT^ . to the lattice data by varying rh and 
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FIG. 5: The eigenvalue density of D5 for /Si^^ = 2.79, V = 20* 
and \v\ = 2 top: mo = —0.184 and bottom: mo = —0.178. 
The smooth curves are the parameter free predictions from 
WCPT (all parameters are fixed from the fit to the v — 0, 
nio = -0.184 data). 



FIG. 6: The dependence of the eigenvalue density of D5 for 
v = on top: the quark mass rh and on bottom: the lattice 
artifacts ag. The central (Red) line in both plots correspond 
to the (Red) analytic curve which is plotted against the data 
in the top panel of Figure [D 



D. Effects of We. and W7 

The effect of Wq < or W7 < is to smear out the 
analytic predictions so that the bumps in the curves be- 
come less pronounced. In Figure [S] we give an example 
of the effect of d| and aj. The more oscillating curve in 
both plots corresponds to the curve in the top panel of 
Figure m which has ag = 07 = 0. As a less pronounced 
structure of the oscillations due to the individual eigen- 
values is unfavored by the lattice data we conclude that 
|dg| and \a^\ appear to be smaller than 0.01. Note that 
Si only enters as squares in the action that is why it 
is natural to give the bound on the square. 



E. Individual eigenvalue distributions 

In order to understand in greater detail whether our 
lattice data match the predictions from WCPT and 
WRMT it is advantageous to consider the distributions of 
individual eigenvalues. As we do not have any simple ex- 
pression available from the chiral Lagrangian approach, 
we will simply make use of the relation to Wilson chi- 
ral random matrix theory, and generate the distributions 
numerically. See Figured 

For the fUjw = 2.635 data on the 16* lattice the first 
eigenvalue matches well that generated from WRMT. 
The second eigenvalue however is somewhat broader. 

In the Piui = 2.79 data from the 20^^ lattice the 
first eigenvalue again matches well that generated from 
WRMT. This time the second eigenvalue is repelled from 
the first in a manner comparable to the one observed in 
the simulation of WRMT. 




FIG. 7: The eigenvalue density of D5 for u — 0. Shown are 
the lattice data and best fit, previously displayed in the top 
panel of Figure [S] along with the a — result at the same 
value of rh. Clearly it is essential to include the order 
effects in WCPT in order to match the data. Note that scale 
on the y-axis is different from that in the other plots. 



F. The real eigenvalues of Dw 

From the flow of as a function of m we have also 
determined the real eigenvalues of Dyy as well as their 
chiralities. This enables us to measure both the density 
of the real modes 

P.^az(A^;«) = ^E'^(^'^-^r)^ (19) 

as well as the distribution 

P^(A^;«) = (E^(^'^ - Ar)sign(fc|75|fc>^ . (20) 

The analytic predictions for both these quantities are 
available, see [3 for preai and [1, Q for p^. The two 
distributions both merge to a ly-fold (5-function at the 
origin for a — > 0. 

From the fit to ps we have already obtained the rel- 
evant values of 'SV, rh and d. This enables a highly 
nontrivial parameter free test of WCPT. In Figure [TT] 
the parameter free analytic predictions are plotted along 
with the results from the lattice. Since accidentally the 
value d = 0.35 was obtained for both data sets we have 
combined the plots. For the plot we have rescaled the lat- 
tice eigenvalues, shown in Figure [TUl by Y,V (the value 
for which we obtained from the respective fits to ps) 
and then shifted the eigenvalues by JlVmo — rh. For the 
20'* data with mo = -0.184 and = 2.79 this reads 
SFmo - rh = 216(-0.184) - 3.5 = -43.2 while for the 
16'' data with /J/^, = 2.635 and mg = —0.216 it reads 
EFmo - m = 157.5(-0.216) - 3.5 = -37.5. 



m=4.S a^-O.O a^=0.0 Oj=0.J5 

m=4.8a=t).0a=i0.2a^^.35 ~ 
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FIG. 8: The eigenvalue density of D5 for v = Q top: the 
dependence on fie and bottom: the dependence on 07. The 
more oscillating (Red) line in both plots has ae = ay = 
and corresponds to the (Red) analytic curve which is plotted 
against the data in the top panel of Figure U Both a% < 
or < tends to smear out the oscillations. Since the 
magnitude of the oscillations in the lattice data presented in 
Figure |4] is of the same magnitude as those in the analytic 
prediction with as = 07 = we can put an upper bound on 
these |ai|, la?] < 0.01. 

We observe that the scaling implied by WCPT is re- 
markable accurate. After the parameter free shift of the 
eigenvalues the center of the data peak is right at the ori- 
gin as predicted by WCPT. Moreover, the width of the 
peak in both data sets is exactly of the right magnitude. 
We also note that the peak height in the data increases 
towards the analytic prediction when going from 16^ to 
20I 

However, in the data there is an asymmetry of preai 
and which persist when going from 16^ to 20^. More- 
over, the average number of real modes is predicted to be 
1.017 by WRMT ior v = I and a = 0.35 0], in the data 
we find 1.237 for the 16'' lattice and 1.185 for the 20^^ 
lattice. That is, the average number of additional real 
modes is off by a factor of 10. We note, however, that 
this average number of additional real modes decrease by 
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FIG. 9: The distribution of the first (black) and second 
(blue) eigenvalue top: in the 16* Piw ~ 2.635 data set and 
bottom: in the 20"* = 2.79 data set. The (red) his- 
tograms show the first and second eigenvalue distribution in 
a WRMT simulation with n = 100 and 100000 matrices. 



22% when going from a = 0.93 fm on 16^ to a = 0.075 
fm on 20"*. 

The asymmetry of the distribution of the real eigen- 
values is a clear signal that at the present volumes there 
are small real eigenvalues of Dw which are not captured 
by leading order WCPT. In [2J| it was observed that the 
eigenvectors corresponding to eigenvalues of Dw which 
fall to the right of the main band have a strong tendency 
to be localized. Obviously such localized modes are not 
described by WCPT. 



VI. CONCLUSIONS 

The analytical predictions for the quenched micro- 
scopic spectral density of the Hermitian Wilson-Dirac 
operator have been compared to the low lying eigenval- 
ues of Wilson fermions on the lattice. In this first study, 
we have demonstrated that the results given for fixed 




FIG. 10: The distribution of the real modes of Dw before 
subtraction of the bare mass and scaling of the axis. 



index of the Wilson-Dirac operator can be matched by 
the predictions from Wilson chiral perturbation theory 
with Ws, > and Wq = Wj = 0. We have successfully 
tested the finite volume scaling at fixed lattice spacing 
as well as the mass scaling as fixed lattice spacing pre- 
dicted by Wilson chiral perturbation theory against the 
lattice data. The spectrum of D5 in this way gives a 
very direct way of measuring the low-energy constants of 
Wilson chiral perturbation theory. Although these con- 
stants represent lattice artifacts, it is essential to know 
their values in order to extract the physical constants. 
Moreover, the detailed predictions can be used with ad- 
vantage to understand how the spectrum can be made to 
retain its gap as simulations are pushed towards physical 
values and towards the continuum limit. 

We emphasize that the assignment of index ly to a, given 
gauge field configuration must be based on the spectral 
flow of the Hermitian Wilson-Dirac in order for the match 
to Wilson chiral perturbation theory to be valid. In the 
microscopic limit the analytic results predict that the 
real eigenvalues of Dw are located in well separated re- 
gions. On the lattice, it should therefore be possible in 
a clean way to introduce a cut-off such that we consider 
only the physical branch of the spectrum. In support of 
this we have explicitly computed the real eigenvalues of 
Dw and found that the bulk of the real eigenvalues fol- 
low the predictions from WCPT. However, we observed 
also an asymmetry of the distribution of the real eigen- 
values not described by WCPT. For this reason some 
configurations that are difficult to classify remain in our 
samples. Similar concerns have been voiced already in 
the context of the overlap operator [2^, where localized 
modes that intuitively should play no role in the contin- 
uum limit nevertheless significantly deform the spectrum 
of the overlap operator at fixed v. We stress that such 
localized modes are not described by Wilson chiral per- 
turbation theory. As such modes are more easily identi- 
fied in the complex part of the spectrum of Dw a careful 
and high-statistics study of the microscopic spectrum of 



10 



Q. 

0.1 











, 1 , 1 1 1 1 1 1 1 1 

Ev_reaL279_20_HYP_QlandQml 






Ev real 2635 16 HYP QlaiidQml 






D=0.J5 






PREDICTION 


- 


I 

1 








V, I^T=^^^&1Q-.u-_l J-_ ^ 1 _ , 1 



10 10 20 30 40 50 60 70 



0.3 




1 1 1 1 1 1 1 1 1 1 1 

Ev_reaL279_20_HYP_Ql 




Ev real 2635 16 HYP Ql 

a=0.35 

PREDICTION 


"S" 0.2 


1 

r 


L 
I 


7f 

>^ 

Q. 






01 









■ \ 








J I \ I I I \ I \ I I I \ I \ I I 

-10 10 20 30 40 50 60 70 



FIG. 11: top: The distribution of the real modes, bottom: 
the distribution defined in Eq. (|20|) . All parameters are 
fixed from the fits to ps . The eigenvalue data are first rescaled 
by TiV and then shifted by moTiV — m where mo is the lattice 
mass and m is the value obtained previously for this mo. 



the Wilson-Dirac operator which address this issue would 
provide a valuable test of Wilson chiral perturbation the- 
ory at presently realistic lattice spacings. The analytic 
prediction for the quenched microscopic eigenvalue den- 
sity of Dw in the complex plane has been presented very 
recently [7|. 

It is also possible to transform all analytical predic- 
tions for the microscopic spectral density of the Hermi- 
tian Wilson-Dirac operator into corresponding expres- 
sions for the pseudoscalar condensate (■ip^^tp), as dis- 
cussed in Ref. [1^. It would be interesting to use this 
quantity to extract the Wilson low-energy constants. 

To summarize we have demonstrated here that the the- 
oretical predictions of Wilson chiral perturbation theory 
with Wq = Wr = and Wg > for the microscopic 
spectral properties of the Wilson-Dirac operator can be 
matched to lattice data. It would be most interesting to 
extend this quenched study to the physical case of two 
light flavors. The detailed theoretical predictions have 
been given recently Q. 

End note: We learned at the Lattice 2011 meeting 
that Albert Deuzeman, Urs Wenger and Jair Wuilloud 
have been investigating the same issues we have discussed 
here. We thank these authors for coordinating publica- 
tion of results. Their paper is Ref. [26l |. 
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